Magnetic friction: Prom Stokes to Coulomb behavior 



Martin P. Maeiera Q Sebast ian Angst, Alfred Hucht, and Dietrich E. Wolf 
Faculty of Physics and CeNIDE, University of Duisburg-Essen, D-47048 Duisburg, Germany 

(Dated: November 11, 2011) 

We demonstrate that in a ferromagnetic substrate, which is continuously driven out of equilibrium 
by a field moving with constant velocity v, at least two types of friction may occur when v goes to 
zero: The substrate may feel a friction force proportional to v (Stokes friction), if the field changes 
on a time scale which is longer than the intrinsic relaxation time. On the other hand, the friction 
force may become independent of v in the opposite case (Coulomb friction). These observations are 
analogous to e.g. solid friction. The effect is demonstrated in both, the Ising (one spin dimension) 
and the Heisenberg model (three spin dimensions) , irrespective which kind of dynamics (Metropolis 
spin-fiip dynamics or Landau-Lifshitz-Gilbert precessional dynamics) is used. For both models the 
limiting case of Coulomb friction can be treated analytically. Furthermore we present an empiric 
expression refiecting the correct Stokes behavior and therefore yielding the correct cross-over velocity 
and dissipation. 
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Friction phenomena, despite their huge importance in 
everyday life, are still not fully understood. Different 
friction mechanisms are possible, leading to different de- 
pendencies of the friction forces on the driving velocity. 
Microscopically, one often assumes Stokes-like friction, 
i.e. a linear velocity dependence. However, this atomistic 
view is in conflict with Coulomb friction at the interface 
between solids, because it approaches a nonzero abso- 
lute value in the limit of small velocities, independent of 
the materials and their surface conditions [THS]. A pos- 
sible solution was offered by the simple model developed 
by Prandtl and Tomlinson, in which a stick-slip insta- 
bility was responsible for Coulomb friction [71 [5] . They 
suggested a surface atom to be coupled by a spring of 
stiffness A: to a slider which moves with constant velocity 
V. The atom interacts with the surface via a periodic 
potential and experiences a viscous friction force propor- 
tional to its velocity x. If k is sufficiently small with 
respect to the potential height, the atom first gets stuck 
in the potential minima and slips when the tension gets 
large enough. The slip motion x does not depend on 
the slider's velocity v, and one observes Coulomb fric- 
tion. However, when k is large with respect to the po- 
tential height, the atom moves with the slider's velocity 
and the friction force is Stokesian. The cross-over from 
one regime to the other has been studied recently |9] . 

What remains a puzzle, however, is that Coulomb fric- 
tion is a far more general phenomenon than one might 
expect from the Prandtl- Tomlinson model, which is for- 
mulated in terms of elastic forces in a periodic potential. 
For example. Stokes as well as Coulomb behavior has also 
been observed for magnetic friction, where elastic forces 
are absent. Being guided by a detailed investigation of 
the crossover between both types of magnetic friction, a 
unifying principle can be formulated that applies to the 
magnetic as well as to the elastic case. 

So far, magnetic friction has been studied in two dif- 
ferent types of models. Ising models with single-spin-flip 



dynamics, where two half spaces move with respect to 
each other, yield Coulomb friction [TUHT2]. Analogous re- 
sults have been obtained in the Potts model [T3 . On the 
other hand, a magnetic dipole scanning a Heisenberg sur- 
face showed Stokesian friction fMUTT] (always provided 
the velocity is not too large). 

Recently a work has been published, in which a point- 
like magnetic perturbation moves through an Ising model 
PEl HH] . The authors claim to have observed Stokes fric- 
tion, which is in conflict with our results for similar mod- 
els [ini [11] . Here we present an explanation of this dis- 
crepancy and clarify, under what conditions Stokes re- 
spectively Coulomb friction occurs. 

The systems studied in Refs. [TUHTB] have in common 
that the motion occurs in a discretized way: The sys- 
tem is at rest for a certain number a/v of Monte Carlo 
sweeps (MCS), after which one half space is moved by one 
lattice constant a. Accordingly we have a periodic exci- 
tation and relaxation procedure, where excitation is fast 
(happens in between two subsequent spin flip attempts) , 
whereas relaxation extends over a/v MCS. By contrast, 
in [TH - IT7] excitation is slow, because due to the dipole- 
dipole interaction a substrate spin feels the approaching 
tip a long time in advance. 

Now we present a simple Id model that interpolates 
between both cases: We consider a position dependent 
field hz{r'), which is moved continuously with constant 
velocity v. r is given in units of a, and v = r. Then 
the discrete motion can be modeled as a step function, 
as shown in Fig. [l] as solid line. For a certain time 1/v 
exactly one spin is exposed to the field with constant am- 
plitude until the field reaches the next spin. Additionally 
the amplitude of the dipole field used in [TiHTT] is plotted. 
From Ref. [TS] we know that for this case the adjustment 
of the spins with respect to the moved field happens in 
an adiabatic way, or in other words the time scale of re- 
laxation is below that of the excitation. To generalize 
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FIG. 1. The dynamics in the different studies can be mapped 
on a time dependent field with amplitude hzir') (here normal- 
ized by its maximal value h) interacting with the spins posi- 
tioned at integer sites. The discrete motion in the Ising/Potts 
model then corresponds to a step-function, which may be 
treated as a fixed spin interacting via exchange with one part- 
ner on the chain. The amplitude of a dipole field is sketched 
for comparison. The field used in this work may be tuned 
by adjusting the parameter 5r from one limiting case to the 
other. 



these setups, we consider a field with steepness 5r <^l, 

h 
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which may be tuned from the step-hke field {Sr = 0, now 
called limiting case (i)) to a slowly varying field {Sr w 0.1, 
case (ii)). By shifting this field according to r' — r — vt 
PP] we can directly influence the time scale at which the 
excitation at a flxed position r occurs, Tgwitch (x 5r/v. 

We first consider a chain of classical, normalized 
Heisenberg spins (|Sr-|=l) of length L with lattice spac- 
ing a, which interact with the field defined above. The 
corresponding time-dependent Hamiltonian is 

L 

m) ^-Y^JSr- Sr+l + d^Sl, + Kir - vt)Sr,., (2) 
r=l 

with the exchange constant J. To get a well defined 
ground state, we use an easy axis anisotropy {dx>Q) and 
anti-periodic boundary conditions S^+l^— S^- The spins 
perform Landau-Lifshitz-Gilbert dynamics [lH |22] , 
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consisting of a precessional motion with a frequency pro- 
portional to "f/ Us, and a damping with the damping con- 
stant a. For simplicity, we neglect temperature here, and 
the dynamic parameters yield a spin relaxation time t^c\- 
The friction force F can be either calculated from the dis- 
sipated power Pdiss or the pumping power Ppump, which 
are equal in the stationary state due to energy conserva- 
tion and therefore we subsequently use F = (P) /v after 
time averaging. The two cases can be described by 
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with the dissipation exponent (f>=l (0=2) for the 
Coulomb (Stokes) case. P can be extracted from the 
energy terms by 



P(i) = 



dn 

dt 
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dhz{r — vt) 
dt 



(5) 



which represents the power pumped into the system by 
the motion. In our simulations (c.f. Fig. [3^) we found 
4>=2 for large Tgwitch, which corresponds to the results in 
P^HTTj . For sufficiently small Tgwitch we get </>=!, which 
was known from simulations in the Ising model, and was 
now reproduced in the Heisenberg model. 

In the following we calculate the velocity Ux at which 
a cross-over from one regime to the other occurs. For 
case (i) only two spins contribute to the sum in Eq. ^ 
at the discrete times vt G Z (at all other times and po- 
sitions the field remains constant), and we can calculate 
the averaged pumping power by discretizing dth^, 
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For the time r^a, corresponding to the time at which the 
amplitude of the field stays nearly constant, no pumping 
(respectively excitation) occurs. We consider rca>Ti.ci, 
i.e. the system always relaxes to equilibrium after a 
pumping event. Since the equilibrium configuration does 
not depend on the dynamics, Eq. |6] tells that here (j)=l. 
The equilibrium configuration for our choice of boundary 
conditions is a domain wall (DW) state, where the out- 
of-axis component is determined by the field and thus 
points in z-direction. As the field interacts mainly with 
only one spin, the shape of the DW is not infiuenced by 
h and we may use the continuum limit profile (a— >0) |23| 
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with the DW width £=\/ J/(2dx), which can be calcu- 
lated from minimizing the free energy [21]. By insert- 
ing (So')=mH(0) and (Si')=m"(l) into Eq. (jej) we now 
can calculate the power which is pumped into the sys- 
tem during each switching event. This quantity can be 
visualized in a potential plot. We again assume that h 
does not infiuence the shape of the DW but its center 
As in limiting case (ii) the system is always near 
equilibrium, and in limiting case (i) it always reaches 
the ground state before being excited out of equilibrium, 
this assumption is justified and we can describe the whole 
configuration with r^^. We look at one cycle at which the 
field's peak moves from to 1, corresponding to the times 
-~l/2<vt<l/2. For given t we can calculate the system's 
total energy as a function of r^w (see the potential lines 
in Fig. [2| . If the system evolved quasi statically, it would 
always be in the current potential minimum. In this pic- 
ture Pc/v corresponds to the energy difference between 
the energy at rdw=0, vt=—l/2 (the equilibrium state) 
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FIG. 2. The system can be parameterized by the center of 
the DW rdw, thus for different times the total energy of the 
system can be calculated. This is done for field parameter 
5r=0.1 for 10 equidistant times (blue to purple curves). If the 
system evolved quasi static, it would follow the configuration 
of minimal energy, marked by the curve. Simulation results 
show both cases, the points indicate energy vs. rdw for 10 time 
steps: (i) ■ The energy Pq /v is periodically pumped into the 
system, which relaxes independent from Tswitch afterwards. As 
the switching occurs at vt — 0, we first see the relaxation from 
the preceding excitation in the left minimum, and at vt = 1/2 
the relaxation in the right minimum is not finished, (ii) • The 
system follows with a lag, but stays near equilibrium, slightly 
above the minimal energy configuration. 



and vt=l/2 (the state which is present when the peak 
of the field has moved to the next spin while the DW is 
still at the same site). Results from simulations (plotted 
as squares in Fig. [2]) confirm this: At vt—0 the system is 
excited to the upper state in a short time, and relaxes to 
the new ground state by adjusting r^jw slowly afterwards, 
until it reaches the new ground state configuration with 
?'dw=l- Simulations of the second limiting case (circles 
in Fig. [2| confirm that the system is always near equilib- 
rium, thus the DW slightly lags behind the ground state. 
From our simulations, we found the pumping power 

dx^s 
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The factor dxv/{5rJ) originates in the synchronization 
with the field, which changes at the time scale Tgwitch- 
The factor av/ig/j emerges from spin dynamics, yielding 
a retardation of the DW as derived in Ref. [U]. Setting 



yields the cross-over velocity 
TT hSr 7 J 

where the system performs a cross-over from the Stokes- 
friction state to the Coulomb-friction state. In Fig. |3] 
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FIG. 3. (Color online) Dissipated power vs. velocity (in natu- 
ral units) for the a) Heisenberg {dx=0.5J,a=0.5,h=J, blue) 
and the b) Ising {(5=1/ J, h=10J, black) model. The simulated 
Sr are Sr=oo (A), 10"'' (T), 10"^ (•), 3 x lO-^(o), IQ-^ (*) 
and 10"^ (□). The grid lines mark the corresponding Vx , the 
dashed lines display the calculated Pc (v) from Eq. ([6| and fit- 
ted Ps{v). For c) we calculated explicitly the cross-over quan- 
tities Px and Vx for both models from Eqs. ( 9|10 1, and plot 
the data again rescaled. Additionally we varied: a=0.3 (pur- 
ple), h=2J (green) and dx=Q.2^J (yellow) for the same 5r set. 
In the inset we plot an effective exponent (j>cfi=d log P/d log v, 
and get a universal cross-over from 2 to 1. 



these cross-over quantities have been calculated and the 
simulation results have been rescaled appropriately. The 
simulation data fit excellent over several magnitudes with 
the derived cross-over quantities, remaining deviations 
are discussed below. We performed also simulations of 
the isotropic Ising model with the same field and peri- 
odic boundary conditions {Sr+N=Sr)- The Ising spins 
undergo spin fiip dynamics with Metropolis probability 
|26) : Randomly chosen spins are fiipped with the prob- 
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ability p/= min (1, exp (— /3A£')), where /3 is the inverse 
temperature and A_E the energy difference between the 
flipped and the not flipped state. 

For (i) we again flnd a behavior ip—1. We assume that 
for (i) the spins relax after each switching event to the 
ground state profile which can be calculated via transfer 
matrix methods [TT] : 



m^(r') = tanh(^/i)[tanh(^J)]l'^'l. 



(11) 



With (S'o',^)=mi(0) and {Sy ^^)=m^{l) in Eq. we get 
=2 is observed for (ii), and we fitted 



Pg oc v^/Sr. 



(12) 



We calculated again the cross-over velocity , which 
is additionally plotted in Fig. [Sjj, and rescaled all data 
points for the cross-over plot, Fig. [3]:. 

Comparing Fig. [3^ and Fig. [Sfj we come to the main 
result of our investigation, namely coincidence concern- 
ing the cross-over between both models, despite the sub- 
stantial remaining differences like e.g. the dynamics of the 
models. The present deviations from the cross-over curve 
are discussed below. The slight increase of P, observed 
in the regime t;>0.1 for all Sr in the Ising model is due 
to the fact that the system has not enough time to relax 
back to equilibrium before the next shift takes place and 
TO^(l) becomes significantly smaller than its equilibrium 
value. As the Heisenberg model contains spin wave exci- 
tations, we observed the generation of spin waves above 
a threshold velocity [TB]. In the cross-over plot these 
spin waves cause a kink above w/fx=0.1 for (and 
a higher peak in the effective exponent plot). For very 
high velocities we observe a lowering of the power, which 
is due to a segregation of the peak of the field and the 
DW, leading to a reduced m^(0)<l. This state with low- 
ered dissipation has already been observed and reported 
in Ref. [H]. 

In conclusion, we presented a new model, which for 
the case of magnetic friction shows a transition from 
Stokes to Coulomb behavior, analogous to the Prandtl- 
Tomlinson model for solid friction. Whereas there, the 
elastic stiffness of the slider was the crucial parameter, 
it is the switching time of the magnetic field in our case. 
The comparison of both models sheds new light on the 
univeral origin of Coulomb behavior, which is based on 
a separation of the relaxation time from the much larger 
time scale on which the system gets excited. Our findings 
are in accordance to field theoretical results by Demery 
et al. [mUnj, who also found Stokes-like friction as their 
system model does not contain discrete sites and thus 
Tea = 0, i.e. the field is continiously driving the system. 
However, their simulation results are not correct as they 
simulated an Ising model with a discontinuous motion of 
a field, which is known to show Coulomb friction. This 
discrepancy stems from an incorrect definition of the fric- 



tion force (Eq. (50) in Ref.|19j. a correct definition has 
been presented in [10]). 
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